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We have measured the thermal conductivity of the one-dimensional (ID) S = 1/2 Heisen- 
berg antiferromagnetic spin system of Sr 2 Cui_ a; Pd a: 03 single crystals including nonmagnetic 
impurities of Pd 2+ . It has been found that the mean free path of spinons along the ID spin 
chain at low temperatures is very close to the average length of finite spin chains between spin 
defects estimated from the magnetic susceptibility measurements. This proves that the ther- 
mal conduction due to spinons at low temperatures in Sr2Cu03 is ballistic as theoretically 
expected [ Zotos et al, Phys. Rev. Lett. 55 (1997) 11029]. 

KEYWORDS: ballistic thermal conduction, thermal conductivity, one-dimensional quantum 
spin system, Sr 2 Cu03 

1. Introduction 

Recently, thermal conductivity in low-dimensional quantum spin systems with the spin 
quantum number S = 1/2 has attracted interest, because the thermal conductivity due to 
spin excitations, K sp m, has been found to be high in various materials, such as the 2deg spin 
ladder system Sr^C^-iCUi, 1-5 one-dimensional (ID) antiferromagnetic (AF) spin systems 
Sr2Cu03, 6 ' 7 SrCu02, 7 Bechgaard salts 8 and BaCu2Si207. 9 It has been found that one of 
essential factors of high K S pin is a large bandwidth of the spin excitations bringing on a high 
velocity of the spin excitations and that the AF correlation between the nearest neighboring 
spins is more suitable for high K sp m than the ferromagnetic one. 10 However, the mechanism 
of K S pi n in low-dimensional quantum spin systems has not been understood fully. 

Accordingly to the theoretical study using the Kubo formula, 11 it has been predicted 
that the thermal conduction due to spin excitations is ballistic at finite temperatures in ID 
spin systems with 5 = 1/2 described by integrable Hamiltonian's, because the heat flow is 
a conserved quantity. 12-16 That is, in these systems, the thermal conductivity due to spin 
excitations, namely, due to spinons possessing 5 = 1/2, K sp i non , is expected to be very high, 
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because the mean free path of spinons, Spinon , is infinite in the ideal case. In spin systems 
described by non-integrable Hamiltonian's, on the other hand, it has been predicted that the 
thermal conduction due to spins is diffusive and that Z sp mon is comparable to the distance 
between the nearest neighboring spins in the high temperature limit. 

As for Sr2Cu03, it has been found from magnetic susceptibility, 17 ' 18 specific heat 18 and 
midinfrared optical absorption measurements 19 that the intrachain exchange interaction, J, 
is as large as more than 2000 K. On the other hand, the interchain exchange interaction, J', 
is as small as ~ 10~ 4 J, according to the estimate from the AF transition temperature, Tn, 
~ 5.4 K 20 ' 21 using the simple relation, Tn ~ \/ J x J' 2 . Therefore, S^CuOs is regarded as an 
almost ideal ID 5 = 1/2 Heisenberg AF spin system described by the integrable Hamiltonian, 



Actually, a high K sp mon has been observed in S^CuOs, as mentioned above. ' The contri- 
bution of Atspinon to the thermal conductivity in S^CuOs is characterized by a clear shoulder 
around 70 K in the temperature dependence of the thermal conductivity along the direction 
parallel to the spin chain only. The existence of the ballistic thermal conduction due to spinons 
has been insisted from the analysis of l sp mon- Besides, it has been reported from the NMR 
measurement of Sr2Cu03 that the spin-diffusion constant related to /-Cspinon is very large, 22 ' 23 
though there is a report that the spin transport is diffusive at finite temperatures. 24 Moreover, 
a similar large spin-diffusion constant has also been obtained from the NMR experiment in 
the ID S = 1/2 AF spin chain ststem a-VO(P0 3 ) 2 . 25 

According to the discussion by Sologubenko et a/. 6 ' their insistence on the existence of 
the ballistic thermal conduction due to spinons in Sr2Cu03 is based on the result that the 
characteristic temperature of the spinon scattering is the order of the Debye temperature, 
0D- It may follow the major scattering is spinon-phonon scattering rather than spinon-spinon 
scattering, but it seems too rough to insist that the thermal conduction due to spinons is 
ballistic because of the absence of the spinon-spinon scattering. Therefore it is not clear 
whether the thermal conduction due to spinons is ballistic or not. In this paper, in order 
to confirm the ballistic nature of the thermal conduction due to spinons in Sr2Cu03, we 
have grown S^Cui-^Pd^Os single crystals with x = 0, 0.004, 0.010 including nonmagnetic 
impurities of Pd 2+ , in which the average length of finite spin chains between spin defects, 
£imp; is expected to decrease with increasing x. We have measured the thermal conductivity 
along the 6-axis parallel to the ID spin chain, k^,, and along the a-axis perpendicular to the ID 
spin chain, K a . We have measured the specific heat to estimate the value of ©d and the specific 
heat of spinons, C sp i non - We have also measured the magnetic susceptibility to estimate the 
amount of spin defects. Then, we have compared the value of Z sp inon estimated from K sp i n0 n and 
Cspinon with L imp estimated from the magnetic susceptibility measurements. As a result, it 
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has been found that the value of i S pinon at low temperatures is very close to the value of L[ mp , 
meaning that spinons are moving along the spin chain between spin defects without being 
scattered at low temperatures. This proves that the thermal conduction due to spinons at low 
temperatures in S^CuOs is ballistic. The preliminary results have already been reported in 

2fi 

our previous paper. 
2. Experimental 

Single crystals of Sr2Cui_ x Pd x 03 with x = 0, 0.004, 0.010 were grown by the Traveling- 
Solvent Floating-Zone (TSFZ) method. In order to prepare the feed rod for the TSFZ 
growth, first, we prepared polycrystalline powder of S^Cui-^Pd^Os by the solid-state re- 
action method. The prescribed amount of SrCOs, CuO and PdO powders with 99.9 % purity 
was mixed, ground, and prefired at 800 °C in air for 24 h. After pulverization, the prefired 
powder was mixed and sintered at 1030 °C in air for 24 h with several times of intermediate 
grinding. After thorough grinding, the powder was isostatically cold pressed at 2.6 kbar into 
a rod of 7 mm in diameter and ~ 120 mm in length. Then, the rod was prefired at 800 °C 
in air for 12 h and sintered at 1040 °C in air for 24 h. As a result, a tightly and densely 
sintered feed rod was prepared. As S^Cui-^Pd^Os melts incongruently, 27 solvent disks with 
the composition of Sr : Cu : Pd = 45 : 55(1 — x) : 55x in the molar ratio were prepared in 
a similar way. The sintering was performed at 800 °C in air for 12 h. The TSFZ growth was 
carried out in flowing O2 gas of 1 bar in an infrared heating furnace. The rotation speed of 
the upper and lower shafts was 15 rpm in the opposite direction. The zone traveling was 1.0 
mm/h. The grown crystals were annealed at 870 °C for 72 h in flowing Ar gas of 1 bar in 
order to remove excess oxygen. The crystals were characterized using x-ray back-Laue pho- 
tography and were confirmed to be a single phase by powder x-ray diffraction. The chemical 
composition of the crystal was determined by inductively coupled plasma optical emission 
spectrometry (ICP-OES) and comfirmed to be the same as the nominal composition. 

Thermal conductivity measurements were carried out by the conventional steady-state 
method. One side of a rectangular single-crystal, whose typical dimensions were 3.5 x 0.8 x 
0.8 mm 3 , was anchored on the copper heat sink with indium solder. A chip-resistance of 
1 kfi (Alpha Electronics Corp. MP1K000) was attached as a heater to the opposite side of 
the single crystal with GE7031 varnish. The temperature difference across the crystal (0.02 
- 1.0 K) was measured with two Cernox thermometers (Lake Shore Cryotronics, Inc. CX- 
1050-SD). The specific heat was measured by the thermal relaxation technique using a PPMS 
(Quantum Design, Model PPMS). The magnetic susceptibility was measured using a SQUID 
magnetometer (Quantum Design, Model MPMS). 
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3. Results and Discussion 

3.1 Magnetic susceptibility 

Figure 1 shows the temperature dependence of the magnetic susceptibility, x, along the 
three principal axes of Sr2Cui_ x Pd x 03 with x = 0, 0.004, 0.010. It is found that \ is almost 
constant at high temperatures and increases with decreasing temperature at low tempera- 
tures below 50 K. The x is so anisotropic that the value of x along the c-axis is larger than 
those along the a- and 6-axes. The increase of x a t l° w temperatures is Curie-like and be- 
comes marked with increasing x. These behaviors of \ are similar to those reported in the 
literature. 17 ' 18,28 

The experimental data of \ are wen fitted using the following equation as in the former 



report: 18 



X = XCuric + Xspin + XO, (2) 



where xcurie is the Curie term, Xspin is the contribution of spin chains and xo is a constant 
term due to the Van Vleck paramagnetism and the ion-core diamagnetism. The xcurie is given 
by 

_ N A g 2 fllS(S +l)x C urie , . 

XCuric " 3k B (T-e) ' [) 

where is Avogadro's number, g the g-factor, /xb the Bohr magneton, xcurie the concentra- 
tion of free spins per Cu, &b the Boltzmann constant and the Weiss temperature. The Xspin 
is given by the following equation proposed by Eggert et al. 29 for the ID S = 1/2 Heisenberg 
AF spin system at low temperatures of k^T J: 

Xspin = ^2J (l + 21n(To/T) ) . ( 4 ) 

where J is the nearest neighbor exchange interaction and To is a parameter depending on 
the second nearest neighbor exchange interaction. In the fitting we assume that To = 7.7 J 
which is true in the case that the second nearest neighbor exchange interaction is neglected. 
The value of the ion-core diamagnetism is put at —10.7 x 10 -5 emu/mol. 30 The value of the 
Van Vleck paramagnetism is fixed at 4.7 x 10~ 5 emu/mol for the a- and 6-axes and 9.7 x 10~ 5 
emu/mol for the c-axis, because its accurate value is not obtained from the fitting. Although 
the fixation of the xo value in the fitting enhances the uncertainty of the obtained value of J, 
it is noted that the change of xo an d J values does not affect the value of the Curie term so 
much, because the Curie term is sensitive to only the increase of the magnetic susceptibility 
at low temperatures. 

The best-fit results are shown by solid lines in Fig. 1. Parameters obtained from the best 
fit are listed in Table I. The value of J is estimated to be 1900 ± 300 K. This value is slightly 
smaller than that in the former report, 18 which may be due to the difference of the fitting 
temperature-range. Values of 6 show no systematic change in correspondence with the Pd 2+ 
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doping, as in the former report. 28 The value of xcurie is found to be about a half of the Pd 
concentration, x. This is reasonably explained as follows. That is, in an ID AF spin system, 
the Curie term originates from free spins due to finite spin segments separated by spin defects 
such as Pd 2+ , lattice defects and unintended impurities. 31 Therefore, there exist finite spin 
segments with a free spin of S = 1/2 composed of an odd number of spins (called odd spin 
segments) and those with S = composed of an even number of spins (called even spin 
segments). An even spin segment is divided by one spin-defect into an even spin segment and 
an odd spin segment. On the other hand, an odd spin segment is divided by one spin-defect 
into two even segments or two odd segments. Therefore, one free spin of S = 1/2 is induced, 
on average, by two spin-defects. 32 Accordingly, it is reasonable that the value of xcurie is about 
a half of x. The average length of finite spin chains, namely, Li mp can be calculated from the 
following equation: 

Limp = ' (5) 



uric 



where xcurie is the average value of xcuric's in each x and c is the lattice parameter along 
the c-axis; c = 3.499 A. Values of L imn are listed in Table III, using the average value of 
^Curie's obtained from the magnetic susceptibility measurements in fields parallel to the three 
principal axes. 
3.2 Specific heat 

Figure 2 shows the temperature dependence of the specific heat, C, of Sr2CuC>3. It is found 
that C exhibits a small peak around 5 K which is in good agreement with Tn observed in the 
neutron scattering 21 and /xSR measurements. 20 ' 21 Since there is no electronic contribution to 
C in Sr2Cu03, C is given by the sum of the magnetic specific heat due to spin chains, namely, 
Cspinon and the specific heat of phonons, C p honon> as the follows: 

C — Cspinon ~\~ Cphonon' (6) 

The C S pi n on is given by the following equation based on the ID AF Heisenberg model with 
S = 1/2 at low temperatures of ksT < J: 33 

C - 2 -^T (7) 

"^spinon — r^j ± j (' J 

where N s is the number of spins. The C p honon is given by the following equation based on the 
Debye model: 

_ 12^ Nk B 3 5 

Cphonon — ggpj 1 + " J , (°J 

where iV is the number of atoms and 5 is the coefficient of the anharmonic term. The data of 
C at low temperatures above Tn are well fitted using Eqs. (6) - (8), as shown by the solid line 
in Fig. 2. Values of the best-fit parameters are J = 2634 ± 150 K, D = 470.8 ± 8.2 K and 
5 = 1.2 ± 0.2 x 10~ 7 J/K 6 mol. These values are comparable with those in the former report, 
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respecively. 6 However, the value of J is a little larger than that estimated from magnetic 
susceptibility measurements, which may be due to the difference of the fitting temperature- 
range bewteen the two kinds of measurement. 
3.3 Thermal conductivity 

Figure 3 shows the temperature dependence of along the spin chain of Sr2Cui_ x Pd x 03 
with (a) x = 0, (b) x = 0.004 and (c) x = 0.010. The temperature dependence of K a perpen- 
dicular to the spin chain in x = is also shown in the inset of Fig. 3 (a). For x = 0, it is found 
that K a increases with decreasing temperature and exhibits a peak around 25 K. This is a 
typical behavior of the thermal conductivity due to phonons, Kphonon- On the other hand, Ky, 
increases with decreasing temperature and exhibits a small shoulder due to /■vspinon 

around 70 

K in addition to the peak due to Kphonon around 25 K. These behaviors are similar to those in 
the former report. 6 ' 7 For x = 0.004, 0.010, both the peak around 25 K and the small shoulder 
around 70 K are suppressed by the Pd 2+ doping. These results indicate that both phonons 
and spinons are scattered by nonmagnetic impurities of Pd 2+ so that their mean free paths 
become short. 

Here, we estimate K S pinon- For this purpose, at first, the estimate of Kphonon is necessary. 
The Kphonon is given by the following equation based on the Debye model. 34 

Kphonon - 27r 2 Uphonon \ h J J Q (<* _ 1)2 ^onon**, W 

where x = tuv/k-QT, uj is the phonon angular frequency, K the Planck constant, w p honon the 
phonon velocity and r p h non the relaxation time of the phonon scattering. The fphonon is 
calculated as 

phonon = ^|£(6vr 2 n)"V3 ) (1Q) 

where n is the number density of atoms. The phonon scattering rate, T~£ onon , is assumed to 
be given by the sum of scattering rates due to various scattering processes as follows, 

ion = ^ + Au,' + Bu'Texp (-^) , (11) 

where Lb, A, B and b are fitting parameters. The first term represents the phonon scattering 
by boundaries; the second, the phonon scattering by point defects; the third, the phonon- 
phonon scattering in the umklapp process. Using Eqs. (9) - (11) and putting 6d at 470.8 
K from the specific heat measurements, the data of K a in x = are well fitted, as shown by 
the solid line in the inset of Fig. 3 (a). The estimated Kphonon in is performed by the fit 
of the data of at low temperatures below 25 K with Eqs. (9) - (11). In the fitting, ©d is 
put at 470.8 K. Values of B and b are put at the same values as these used for the fit of K a 
in x = 0, respectively, because the phonon-phonon scattering in the umklapp process seems 
neither to be affected by the direction nor by the slight doping of Pd 2+ so much. The adjusting 
parameters are only and A, which depend on the phonon scattering by boundaries and 
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by point defects, respectively. Then, K S pinon is estimated by subtracting the fitting curve of 
Kphonon from the data of k D ; as shown by dashed lines in Fig. 3. Here, it is noted that /c S pinon is 
a little underestimated, because K sp i non is neglected at low temperatures below 25 K. Values 
of the parameters used for the best fit are listed in Table II. These are comparable with those 
in the former report. 3 ' 6 ' 7 ' 9 It is found that the value of A increases with increasing x. This is 
resonable, because phonons are scattered by substituted Pd 2+ ions. 

Next, we Gst Illicit 6 ^spinon 

using the following equation, 

^spinon = Cspinon ^spinon ^spinon j (12) 

where w sp inon is the velocitiy of spinons. The f sp inon is given by the following equation based 
on the des Cloizeaux-Pearson mode at low temperatures of k-QT <C J: 35 

^spinon — ^ , {^J 

where a is the distance between the nearest neighboring spins in the chain. Therefore, Spinon 
is calculated using Eq. (7) as follows, 

Zspinon = vriV s afc|T Kspinon - (14) 
Figure 4 shows the temperature dependence of Spinon obtained thus and L[ mp estimated from 
the magnetic susceptibility measurements. It is found that Z sp i non increases with decreasing 
temperature and seems to be saturated at low temperatures. Since K sp i non is neglected at 
low temperatures below 25 K in this analysis as mentioned above, values of Z S pinon at low 
temperatures are uncertain. Therefore, values of l S pinon at low temperatures are estimated as 
shown by solid lines in Fig. 4, by fitting the data of Spinon at high temperatures above 55 K 
with the following simple equation: 6, 7,9 

Unon = ^Texp(-T7T)+L- 1 , (15) 

where A s , T* and L are fitting parameters. The first term is due to the spinon-phonon and/or 
spinon-spinon scattering in the umklapp process and T* is the characteristic temperature. The 
second term is due to the spinon scattering by spin defects. The value of L correspounds to 
the saturated value of Z sp inon at low temperatures. Parameters obtained from the best fit are 
listed in Table III. It is found that the value of ^4 S inreases with increasing x. Since the Pd 2+ 
doping is guessed to induce local phonons scattering spinons around Pd 2+ , the spinon-phonon 
interaction may increase with increasing x. To our surprise in Table III, values of L decrease 
with increasing x and are very close to those of L imp , respectively. This means that / sp i n0 n 
at low temperatures is approximately limited by L imp . That is, the thermal conduction due 
to spinons is limited at low temperatures only by scattering by spin defects. Accordingly, it 
is concluded that the thermal conduction due to spinons at low temperatures is ballistic as 
theoretically expected. 12-16 
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4. Conclusion 

We have measured the thermal conductivity, magnetic susceptibility and specific heat of 
Sr2Cui_ x Pd x 03 single crystals with x = 0, 0.004 and 0.010 which are regarded as an ID 
S = 1/2 Heisenberg AF spin system described by an integrable Hamiltonian, in order to 
prove the theoretical prediction that the thermal conduction due to spinons is ballistic. In 
this system, the length of the spin chain between spin defects is controlled by the doping of 
nonmagnetic impurities of Pd 2+ . We have estimated the average length of finite spin chains 
between spin defects, L[ mp , from the magnetic susceptibility measurements and the mean 
free path of spinons, / sp mon, from the thermal conductivity measurements using the Debye 
temperature estimated from the specific heat measurement. It has been found that values 
of 4pinon at low temperatures for x = 0, 0.004 and 0.010 are very close to those of Li mp , 
respectively. This means that spinons carry heat along the spin chain between spin defects 
without being scattered at low temperatures. Accordingly, our results strongly support the 
theoretical prediction that the thermal conduction due to spinons in ID S =1/2 spin systems 
described by intgrable Hamiltonian's is ballistic. 
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T (K) 



Fig. 1. Temperature dependence of the magnetic susceptibility, \> °f S^Cui-zPd^Os with x = 0, 
0.004, 0.010 in a magnetic field of 1 T parallel to (a) the a-axis, (b) 6-axis and (c) c-axis. Solid 
lines indicate the best-fit results using Eqs. (2) - (4). 
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Fig. 3. Temperature dependence of the thermal conductivity along the &-axis parallel to the spin 
chain, k^, of Sr 2 Cui_ :E Pd 2 ;03 with (a) x = 0, (b) 0.004, (c) 0.010. The inset shows the temperature 
dependence of the thermal conductivity along the a-axis perpendicular to the spin chain, n a , of 
Sr 2 Cu03. Solid lines are Kphonon estimated using Eqs. (9) - (11) on the Debye model. Dashed lines 
are K sp inon obtained by subtracting Kphonon from Kb- 
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T (K) 

Fig. 4. Temperature dependence of the mean free path of spinons, Z sp inon, of S^Cui-zPd^Os with 
x = 0, 0.004, 0.010. Solid lines are the best-fit results using Eq. (15). Dashed lines are the average 
length between spin defects, L[ mp , estimated from the magnetic susceptibility measurements for 
x = (upper), 0.004 (middle), 0.010 (lower). 
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Table I. Parameters used for the fit of the temperature dependence of the magnetic susceptibility, \, 
in SraCui-^Pd^Oa with Eqs. (2) - (4). 



X 


axis 




0(K) 


J(K) 




a 


0.0013(1) 


-3.54(14) 


1680(7) 





b 


0.0011(1) 


-2.99(4) 


1780(2) 




c 


0.0007(1) 


-1.97(13) 


1800(7) 




a 


0.0024(1) 


-0.35(12) 


1650(44) 


0.004 


b 


0.0028(1) 


-0.71(9) 


1720(29) 




c 


0.0027(1) 


-0.65(7) 


2010(37) 




a 


0.0040(1) 


-0.88(6) 


1980(40) 


0.010 


b 


0.0039(1) 


-0.74(6) 


2230(51) 




c 


0.0046(1) 


-0.77(6) 


2220(62) 



Table II. Parameters used for the fit of the temperature dependence of the thermal conductivity, k, 
in SraCu^Pd^Oa with Eqs. (9) - (11). 



X 


axis 


Lb (m) 




A( 


s 3 ) 




B (s/K) 


b 





a 


6.00 x 10" 


-4 


1.38 x 


io- 


-44 


2.70 x 10~ 18 


3.49 




b 


6.00 x 10" 


-4 


1.38 x 


io- 


-44 


2.70 x 10~ 18 


3.49 


0.004 


b 


9.00 x 10" 


-4 


3.96 x 


10" 


-44 


2.70 x 10" 18 


3.49 


0.010 


b 


5.90 x 10" 


-4 


4.57 x 10" 


-44 


2.70 x 10" 18 


3.49 



Table III. Parameters used for the fitting of the temperature dependence of the mean free path 
of spinons, Z sp i non , with Eq. (15) in S^Cui^Pd^Oa. The average value of the concentration of 
free spins per Cu, rrcurie, an d the average length between spin defects, Li mp , estimated from the 
magnetic susceptibility measurements are also listed. 



X 


2-Curic 


Amp (A) 


L(A) 


(s/K) 


T* (K) 





0.0010 


1960 


2110 ±26 


8.89 ±0.47 x 10~ 5 


227 ±5 


0.004 


0.0026 


753 


821 ± 83 


10.1 ± 1.8 x 10~ 5 


136 ± 22 


0.010 


0.0042 


466 


498 ± 60 


20.9 ±4.2 x 10" 5 


145 ± 24 
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